The chromo-Weibel instability in an expanding background 
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In this proceedings contribution we review recent calculations of the dy- 
namics of the chromo-Weibel instability in the quark gluon plasma. This 
instability is present in gauge theories with a one-particle distribution func- 
tion which is momentum-space anisotropic in the local rest frame. The 
conditions necessary for triggering this instability can be present already 
in the color-glass-condensate initial state or dynamically generated by the 
rapid longitudinal expansion of the matter created in a heavy-ion collision. 
Using the hard-loop framework we study the case that the one-particle dis- 
tribution function possesses an arbitrary initial momentum anisotropy that 
increases in time due to longitudinal free streaming. The resulting three- 
dimensional dynamical equations for the chromofield evolution are solved 
numerically. We find that there is regeneration of the longitudinal pressure 
due to unstable plasma modes; nevertheless, the system seems to maintain 
a high-degree of momentum-space anisotropy. Despite this anisotropy, we 
find that there is rapid longitudinal thermalization of the plasma due to 
the non-linear mode couplings inherent in the unstable evolution. 
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One outstanding question in the theoretical study of ultrarelativistic 
heavy ion collisions is the timescale for and processes involved in the ther- 
malization and isotropization of the quark-gluon plasma (QGP). Empirical 
evidence in favor of fast thermalization and isotropization of the QGP gener- 
ated in heavy ion collisions was provided by the success of phenomenological 
relativistic hydrodynamical models [l}|9] . The success of these models in de- 
scribing the collective flow observed at the Relativistic Heavy Ion Collider 
(RHIC) and the Large Hadron Collider (LHC) suggests that the QGP may 
become thermal and isotropic on rather short time scales. However, in recent 
years there has been an important realization that successful phenomeno- 
logical application of viscous hydrodynamics may not necessarily imply fast 



isotropization of the QGP in heavy ion collisions [5 10-15 . Currently, the 



question of the degree of momentum-space isotropy of the QGP generated 
in heavy ion collisions is an open question. In this paper we review recent 
numerical calculations |16| which utilize the hard-thermal-loop framework 
description of an anisotropic QGP. 

Due to the rapid longitudinal expansion of the quark gluon plasma, one 
expects generation of momentum-space anisotropics in the pt-Pl plane. In 
the weak-coupling limit the system is expected to be highly-anisotropic at 
early times. In weakly-coupled quantum chromo dynamics (QCD) the pres- 
ence of momentum-space anisotropics induces unstable plasma modes. The 
existence and properties of these unstable modes has been studied using ki- 



netic theory and diagrammatic methods 17-27 . This instability has been 



dubbed the chromo- Weibel instability in reference to the analogous Weibel 



instability which exists in Abelian electromagnetic plasmas 28 . In the 
weak-field regime with a fixed momentum-space anisotropy, the chromo- 
Weibel instability initially causes exponential growth of transverse chromo- 
magnetic and chromoelectric fields; however, due to non-Abelian interac- 
tion between the fields, exponentially growing longitudinal chromomagnetic 
and chromoelectric fields are induced which grow at twice the rate of the 
transverse field configurations. Eventually, all components of the unsta- 
ble gauge- field configurations become of equal magnitude. As a result, one 
finds strong gauge field self-interaction at late times and numerical simu- 
lations are necessary in order to have a firm quantitative understanding of 



the late-time behavior of the system 23 , 29 - 44 



In order to understand the precise role played by the chromo- Weibel 
instability in ultrarelativistic heavy ion collisions one must include the ef- 
fect of the strong longitudinal expansion of the matter. For the first few 
fm/c of the QGP's lifetime the longitudinal expansion dominates the trans- 
verse expansion. Therefore, to good approximation, one can understand 
the early time dynamics of the quark gluon plasma by considering only lon- 
gitudinal dynamics. The first study to look at the effect of longitudinal 
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expansion was done in the context of pure Yang-Mills dynamics initialized 
with color-glass-condensate initial conditions onto which small-amplitude 
rapidity fluctuations were added (34j. The initial small-amplitude fluctua- 



tions result from quantum corrections to the classical dynamics [35 43 , 45] . 
Numerical studies have shown that adding spatial-rapidity fluctuations re- 
sults in growth of chromomagnetic and chromoelectric fields with amplitudes 
~ exp(2m£) y/r/Qs) where mP D is the initial Debye screening mass and r is 
the proper time. This growth with exp(y'r) was predicted by Arnold et al. 



based on the fact that longitudinal expansion dilutes the density 23 



In this proceedings contribution we briefly review our recent paper |16| in 



which we utilized the hard-expanding-loop framework [46 47 to numerically 
determine the evolution of the chromoelectric and chromomagnetic fields 
induced by fluctuations of a system of high-momentum particles which are 
undergoing longitudinal free streaming. Due to the fact that the hard parti- 
cles are longitudinally free streaming their local rest frame momentum-space 
anisotropy £ = ^(Pt) I (Pl) ~ 1 increases as £ = (t/t 1so ) 2 — I, where t 1so is the 
proper time at which the distribution function is assumed to be isotropic. 
In Fig. [T] we plot the unstable 
mode growth rate T/mr> for fixed 
£ as a function of k z /mo where 
m_o is the Debye mass at the 
proper time t; so . As can be seen 
from this figure, as the degree 
of momentum-space anisotropy in- 
creases, more and more modes 
become unstable. Therefore, 
when one has a momentum-space 
anisotropy which is increasing in 
time, more and more modes be- 
come unstable as time progresses. 
In fact, one finds that at late times 
~ rriD^r/Tiso 
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Fig. 1. Unstable mode growth rate T/tud 
for fixed £ as a function of k z /mr) where 
m d is the Debye mass at the proper time 



where k z>raax 
is the wavenumber of the highest 
mode which is unstable. However, 
because of the dilution of the parti- 
cle density due to the longitudinal 

free streaming one finds that the maximal unstable growth rate decreases 
with time as T* ~ mDyThopr. These two effects compete with one another, 
with the former causing unstable growth at higher and higher wave numbers 
as time progresses and the later causing the late time growth to change from 
a pure exponential to exp (2mD^TT mo ) . Both effects should be taken into 
account by the dynamical framework for chromofield evolution. 
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1. Dynamical Framework 

We assume that the background particles are longitudinally free stream- 
ing and, as a result, the background (hard) particles possess a local rest 
frame momentum-space anisotropy which increases monotonically in proper- 
time as specified above. Given an isotropic distribution /i SO , the correspond- 
ing free-streaming distribution is fo(p,x) = fiso{\Jp\ + (p' z t/t iso ) 2 ^ 

Ziso^yVi +P^/ 7 ilo)- Following 46 we obtain the dynamical equation 



obeyed by color perturbations Sf a of a color- neutral longitudinally free- 
streaming momenta distribution /o which can be written compactly as 
V • D 5f a \ pfJ , = gV^F^ 1/ dY p ^fo(p±,Pr])- This equation must be solved simulta- 
neously with the non-Abelian Yang-Mills equations which couple the color- 
charge fluctuations back to the gauge fields via the induced color-currents 



J a 



D,Fr=f a =gt R J J^3^/ a (p,x,t), (1.1) 

where D a = d a — ig[A a ,-] is the gauge covariant derivative and F a p = 
d a Aj3 — dpA a — ig[A a , Ap) is the field strength tensor, and g is the strong 
coupling. The above equations are then transformed to comoving coordi- 
nates with the metric ds 2 = dr 2 — dx^_ — T 2 dr] 2 . 

The resulting dynamical equations are numerically solved in tempo- 
ral axial gauge on a spatial lattice. In order to maintain gauge invari- 
ance with respect to three-dimensional gauge transformations, the spatially- 
discretized fields are represented by plaquette variables and evolved along 
with the conjugate momentum using a leap-frog algorithm. The fluctuation- 
induced currents are represented by auxiliary fields which are discretized in 
space and also on a cylindrical velocity-surface spanned by azimuthal veloc- 
ity and rapidity. As a result, the simulations are effectively five-dimensional 
and are therefore computationally intensive. For details concerning the nu- 



merical implementation we refer the reader to Ref. 16 . 



2. Results 

We used a five-dimensional lattice size of (iVj. x A^) x (N u x N^) = 
(40 2 x 128) x (128 x 32) with transverse spatial lattice spacing of a = Q~ l 
and longitudinal spatial lattice spacing of a v = 0.025. Here Q s is the nuclear 
saturation scale which is approximately 2 GeV and 1.4 GeV at LHC and 
RHIC energies, respectively. For the initial conditions we seeded current 



fluctuations of amplitude A which had a UV spectral cutoff (see Ref. 16 
for details of the spectrum of initial fluctuations). In Fig. [2] (left) we show 
the various components of the chromofield energy density as a function of 
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T = Q s T/10 T 



Fig. 2. On the left we plot the various components of the chromofield energy density 
as a function of proper time. On the right we plot of the total (field plus particle) 
longitudinal over transverse pressure as a function of proper time. 

rescaled proper time f. For LHC and RHIC initial energy densities one unit 
in f corresponds to approximately 1 fm/c and 1.4 fm/c, respectively. For 
this figure an initial fluctuation amplitude of A = 0.8 was chosen. As can 
be seen from this figure after approximately 1 fm/c we begin to see rapid 
growth of the transverse chromomagetic field, followed by the transverse 
chromoelectric field, and then the longitudinal chromofields. In Fig. [2] (right) 
we show the resulting ratio of the total (particle plus field) longitudinal 
pressure divided by the total transverse pressure for various values of A. At 
early times, prior to unstable mode growth, one observes from this figure 
that the longitudinal pressure drops due to the longitudinal free streaming 
of the hard particle background; however, when the unstable modes begin 
to grow, one observes a regeneration of the longitudinal pressure by the 
unstable modes which have their wave vectors pointed primarily along the 
longitudinal direction. In addition, one observes that the time at which 
isotropy is restored is primarily sensitive to the initial fluctuation amplitude 
A. 

In addition to extracting information about the energy density and pres- 
sures of the system as a function of proper time, one can also extract in- 
formation from the gauge field spectra. The longitudinal spectra can be 
obtained following Ref. [45] by first Fourier transforming each field compo- 
nent E±(x±,r]), E r) (x±,r)), B±(x±,rj) and B v (x±,r]), integrating over the 
transverse wave vectors and decomposing each according to the longitudinal 
wave vector u, in terms of which the electric and magnetic energy densities 
are decomposed into longitudinal energy spectra (see Ref. [16] for details). 
One problem with such spectra is that they are not gauge invariant. As 
an additional spectral measure we also extract the transverse momentum- 
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Fig. 3. On the left we plot the longitudinal spectra at various proper times. On 
the right we plot the extracted longitudinal temperature which was obtained by a 
fit (see text) to the longitudinal spectra (£ ) or the Fourier-transform of the spatial 
energy density (£). 

averaged longitudinal spectra obtained by Fourier-transforming the spatial 
distribution of the total field energy density. In Fig. [3] (left) we show the 
extracted longitudinal spectra extracted using the first method averaged 
over 50 runs. The spectra extracted using the second method have similar 
features to the left panel but, due to limited space, we do not show them 
here (see the left panel of Fig. 4 in Ref. [48] for this plot). In Fig. [3] (right) 
we plot the gauge-field temperature extracted from the spectra via fits to 
the form £ tx / dk z (k 2 z + 2\k z \T + 2T 2 ) exp {—\k z \/T) which is obtained by 
integrating a Boltzmann distribution over transverse momenta. In the fig- 
ure we show the fitted temperature obtained from both types of extracted 
spectra (the first method is indicated as '2V and the second method as 
1 Tl (£)')• I n both cases one sees that after an initial period of cooling, the 
gauge sector begins to heat up with the temperatures extracted using the 
two methods being approximately the same. We note that the quality of the 
fits is extremely good (see Fig. 10 of Ref. [16] for comparisons of the fitted 
function to the data at various proper times). The fit function above begins 
to describe the observed spectra very well at early times corresponding to 
f ~ 1 indicating extremely fast longitudinal thermalization of the spectra 
even though the system is still highly anisotropic at this moment in time. 

3. Conclusions 

In this proceedings contribution we have briefly reviewed the recent find- 
ings of our three-dimensional hard-expanding-loop simulations. The chief 
results were: (i) one sees regeneration of the longitudinal pressure by unsta- 
ble chromofield modes, however, the system remains anisotropic for many 
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fm/c; (ii) despite being anisotropic, there appears to be a rapid longitudi- 
nal thermalization due to non-linear mode couplings induced by unstable 
mode growth. In the future we are planning to improve our numerical re- 
sults by utilizing much larger lattice sizes and also studying pure Yang-Mills 
dynamics in an expanding metric. 
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